Phase-matched coherent hard x-rays from relativistic high-order harmonic generation 
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High-order harmonic generation (HHG) with relativistically strong laser pulses is considered em- 
ploying electron ionization-recollisions from multiply charged ions in counterpropagating, linearly 
polarized attosecond pulse trains. The propagation of the harmonics through the medium and the 
scaling of HHG into the multi-kilo-electronvolt regime are investigated. We show that the phase 
mismatch caused by the free electron background can be compensated by an additional phase of the 
emitted harmonics specific to the considered setup which depends on the delay time between the 
pulse trains. This renders feasible the phase-matched emission of harmonics with photon energies 
of several tens of kilo-electronvolt from an underdense plasma. 
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In the last decades, atomic high-order harmonic gen- 
eration (HHG) in the non-relativistic regime pQ has been 
developed to a reliable source of coherent extreme ultra- 
violet (XUV) radiation and of attosecond pulses opening 
the door for attosecond time-resolved spectroscopy 0. 
The further advancement of this technique into the hard 
x-ray domain would, in particular, allow for ultrafast co- 
herent diffraction imaging of single particles, clusters and 
biomolccules with sub-angstrom resolution, tracking the 
electron motion in atoms and even for the investigation of 
time-resolved dynamics of nuclear excitations. The large 
scale x-ray free electron lasers are likely to fulfill this task 
partly but are limited to energies around 10 keV. 

Is it possible to extend the table-top HHG sources into 
the hard x-ray domain? In principle, the harmonic pho- 
ton energy can be increased by using stronger driving 
laser fields. The state-of-the-art technique allows to gen- 
erate coherent x-ray photons up to the keV energy range 
[3] and to produce short XUV pulses of less than 100 as 
[4] from non-relativistic HHG in an atomic gas medium. 
However, progress in this field appears to have reached a 
limit. Most importantly, the further increase of the driv- 
ing field intensity transfers the interaction regime into the 
relativistic domain where the drift motion of the ionized 
electron in the laser field propagation direction prohibits 
the recollision and, consequently, suppresses HHG [5]. 
This happens at laser intensities above of 4x 10 16 W/cm 2 
at infrared wavelengths, corresponding to the HHG cut- 
off frequency of uj c w 10 keV. This indicates the limit of 
non-relativistic HHG. The second point hindering HHG 
at high intensities is the phase-matching problem. In 
strong laser fields, outer-shell electrons are rapidly ion- 
ized and produce a large free electron background causing 
phase mismatch between the driving laser wave and the 
emitted x-rays. 

Various methods for counteracting the relativistic drift 
have been proposed. However, no universal solution has 
been found, each method has its drawbacks. To sup- 
press the drift, highly charged ions moving relativisti- 
cally against the laser propagation direction [B] or a gas 
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FIG. 1: (Color online) The HHG setup with counterpropagat- 
ing APTs: (a) The trajectory of the rescattering electron; (b) 
The contribution of different parts of the medium to HHG. 
The harmonics emitted from the green (dark) area propagates 
along the red (wavy) arrow. The HHG from shaded area is 
damped. The two driving APTs are shown in different blue 
color shades. The arrows indicate the propagation direction. 



of positronium atoms [7] can be used. Different combi- 
nations of laser fields also have been proposed for this 
purpose such as a tightly focused laser beam [S], two 
crossed laser beams with linear polarization [9J [10] or 
with equal-handed circular polarization |llj . In the lat- 
ter field configuration, the relativistic drift is eliminated, 
however, in this scheme the phase-matching is especially 
problematic to realize [T2]. We have shown in [T51 - fT5] 
that strong attosecond pulse trains (APTs) employed as 
a driving field for HHG can be very useful to suppress 
the relativistic drift. However, all these efforts have only 
addressed HHG of a single atom rather than coherent 
emission from a macroscopic gas target. 

In this Letter, we investigate the feasibility of phase- 
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matched harmonic emission from an underdense plasma 
of multiply charged ions for a relativistic HHG setup em- 
ploying two counterpropagating APTs. We show that 
the HHG driven by counterpropagating APTs has an ad- 
ditional intrinsic phase depending on the time delay be- 
tween the pulses as well as on the pulse intensity. This 
phase avails to compensate the phase mismatch between 
the driving laser field and the emitted harmonics due to 
the free electron background. The latter can be achieved 
by modulating the driving field intensity with a slowly 
decreasing envelope. We have performed a complete, 
quantitative analysis of the macroscopic yield of the rel- 
ativistic HHG evidencing a small but detectable signal. 

The applied setup for relativistic HHG is shown in 
Fig. [lja). The driving fields are two counterpropagating 
APTs consisting of 100 as pulses with a peak intensity of 
the order of 10 19 W/cm 2 and a spectral range of about 
20eV. Such pulses can be generated by employing the 
relativistic oscillating mirror of an overdense plasma sur- 
face in a strong laser field [IB]. The electron is liberated 
by one pulse, it is driven in the continuum and under- 
goes the relativistic drift. A moment later, the second 
pulse reaches the electron, it reverts the drift and real- 
izes rescattering. The drift compensation is very efficient 
as one can deduce from Fig. [2] where the single-atom 
spectral emission rate in the present setup is compared 
with that in the dipole approximation (DA). The setup 
exhibits no significant drift any more. The rate for the 
applied setup is much higher than the one for a con- 
ventional laser field with the same cutoff. In the latter 
case, due to the drift, the rate would drop rapidly with 
increasing laser intensity. 

We continue to estimate the macroscopic HHG yield 
of the given setup. The contributions of different parts 
of the medium to the harmonic emission are shown in 
Fig. [ljb) . Emission into the right direction is considered. 
HHG occurs only in small zones of the medium shown as 
dark areas in Fig. |TJb) . In this case, necessarily, re- 
combination has to be arranged by a pulse propagating 
to the right. For instance, ions between pulse 1 and 3 
in the sketch could have been ionized by pulse 1, which 
triggered the process. In the next step, pulse 2 could ar- 
range rescattering and HHG. This is possible only up to 
a certain point (dark area) since then pulse 2 and 3 meet 
and the process starts again. Contributions of atoms ex- 
periencing two pulses simultaneously (light shaded areas) 
are frustrated due to the chaotic trajectories of ionized 
electrons in this region j!7j . This limits the volume in 
longitudinal direction to about 1/3 and the possible de- 
lays are within 1 — 2 fs. After half a period, 1/3 of the 
white area would emit into the right direction. 

In the following, we specify our model. As HHG 
medium, an underdense plasma of 6+ ions (ionization 
potential I p = 27.18 a.u.) is used which is immediately 
formed when the first laser pulse of relativistic intensity 
is applied to a neutral atomic gas. This is because the 



outer shell electrons of an oxygen atom are almost in- 
stantaneously ionized due to a much smaller binding po- 
tential (0.5 a.u. - 5.1 a.u.) other than the two remain- 
ing electrons in the closed ls-shell. HHG is produced 
only by the tightly-bound inner electron having an ioniza- 
tion potential corresponding to the tunneling condition 
at relativistic intensities. The 6+ emission is slightly re- 
duced by the depletion to 7+ whereas the 7+ emission 
is not phase-matched in the proposed phase-matching 
scheme. The driving laser pulses are plane waves numer- 
ically propagated in the relativistic free electron back- 
ground using a Crank-Nicolson-algorithm. The density 
of the free electrons is assumed to be constant because 
the outer shell ionization time is small compared to the 
laser period. Absorption of the high-frequency HHG pho- 
tons can be neglected because their energy is much higher 
than the largest atomic transition energy. In order to find 
the overall HHG yield, the photon spectral density Jrj^- 
in the far field is calculated (in atomic units) [18] : 
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where n' is the emission direction, R the radius at an 
observation point, ujjj the positive harmonic frequency. 
The spectral component of the electric field reads 



i " Rc2 / d 3 x a j a (x a ,u) H ,n'), (2) 



where p is the density of the uniformly distributed ions, 
j a (x,t,x a ) the current density at a space-time point x 
and t of a single ion located at x a and j a {x ai uj H ,n') 




FIG. 2: (Color online) Single-atom HHG rates via [15] in laser 
propagation direction: (red, largest cutoff) for the discussed 
setup; the delay of pulses is 1.5 fs, the laser field strength 
E = 21 a.u., and I p = 27.18 (0 6+ ); (blue, lowest cutoff) 
for the same setup and parameters in DA including the mass 
shift; (black, bottom) for a conventional propagating laser 
field with Eo = 2.7 a.u., I p — 7.35 a.u. (hydrogen-like ion), 
and (gray, top) the latter within the DA including the mass 
shift. The indicated parameters are chosen such that both the 
cutoffs and the average ADK-tunneling rates are the same for 
the two fully-relativistic curves. 



3 



its Fourier transform which is calculated quantum- 
mechanically using the relativistic strong field approxi- 
mation (SFA) [19j . The emitted harmonics have the same 
polarization direction (along the x-axis) as the incident 
laser field. The x-componcnt of the spectral electron cur- 
rent density at a single multiply charged ion in the SFA 
based on the Klein- Gordon equation reads |20) : 

j a (x a ,w H ,n') = -^J^: j ^ x J d A x'(j)*(x ~ x a ,t) 

xV H (x)G(x,x')K{t')V AI (x')<P(x' - x a ,t'), (3) 

where (j>(x — x a ) is the wave function of the bound elec- 
tron, Vh(x) the interaction Hamiltonian of the electron 
with the harmonic field, Vai(x) the potential of the ionic 
core and G(x, x') is the Green function describing the free 
electron evolution in the counterpropagating laser fields 
defined in [TS] . We included a tunneling correction factor 
n(t) = \J wadk (t) Jwk (t) upgrading the Keldysh tunnel- 
ing rate Wif(i) comprised in the SFA to fit the ADK- 
ionization rate [21] , As discussed, we only consider the 
relevant scenario of interaction where the electron moves 
in different counterpropagating pulses successively. Ac- 
cordingly, in each stage of the excursion, we approximate 
the Green function by the Volkov Green function in a 
field of the appropriate single laser pulse (see [H]). Eq. 
([3]) is evaluated in the saddle point approximation. 

Let us have a closer look at the phase difference of 
the harmonics emitted from different ions separated by a 
distance Az in the propagation direction: 

Atp = A arg Ja » Az y - - — (4) 

\ c v g dz ) 

The first term describes the phase mismatch due to the 
free electron dispersion with Av g = oj^/2ojJj, lu p being 
the electron plasma frequency and v g the group veloc- 
ity of the driving laser pulse, whereas the last term is 
the single-atom emission phase (<pi) depending on the 
laser field conditions. This intrinsic phase is deter- 
mined by the classical action of the electron trajectory 
recolliding with the specific harmonic energy and can be 
estimated as (fi ~ U p (r a )T(r a ), with an effective pon- 
deromotive potential U p {r a ) and electron excursion time 
r(r a ). Thus, (/?, depends on the laser intensity as well as 
on the delay between the two pulses. The latter, being 
unique for this laser setup, mainly affects the electron 
excursion time r(r a ) and varies along the propagation 
direction. In order to achieve phase-matching, one can 
vary the laser intensity along the propagation direction 
to balance the intrinsic phase with the phase slip due 
to dispersion. The required intensity variation to have 
a constant complex phase argj a in the entire medium 
is calculated numerically and shown in Fig. [3] for the 
first interaction zone. It is optimized for the long tra- 
jectory of 50 keV energy but could be accomplished for 
any energy value below. Note that only one of the short 
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FIG. 3: (a) shows the optimized peak electric field variation 
in the first interaction zone. The excursion time increases 
to the right, (b) Emitted spectral photon number (spikes at 
integer multiples of uj do not appear due to local averaging) 
of the setup for the medium described in the text. 



and long trajectories can be phase-matched since their 
classical actions are different. For the analytical descrip- 
tion of the spatial variations of the laser field in the ex- 
pression for G(x, x'), the eikonal- Volkov approximation 
is applied [22] . This is justified because the additional 
driving field causing the modulation perturbs the elec- 
tron energy only slightly. The experimental realization of 
the phase-matched scheme could be achieved, e.g., with 
a modulated hollow core waveguide. 

We employ a medium length as short as the spatial 
extent of the APT to minimize dispersion. In our sim- 
ulation, each APT consists of 15 pulses with an APT 
duration of 40 fs. Our calculations show that in the case 
of longer APTs, the pulses in the train strongly spread 
due to dispersion and overlap, thus, violating the con- 
dition for the drift compensation. All pairs of pulses 
have almost the same coherent contribution to the over- 
all yield. Since the pulses in different zones have experi- 
enced a different propagation length through the plasma, 
their shapes differ slightly. However, phase-matching still 
can be maintained by slightly adjusting the modulation 
profile, as long as the pulse shape still supports the rec- 
ollision scheme. The phase-matching scheme imposes a 
strong demand on the jitter of the laser field AE/E: 
AUpT < 1 yields AE/E < ([/p-r)" 1 ~ 10" 4 . We choose 
a gas density of p — 10 19 /cm 3 (ionized by the laser as 
described before), a diameter of 1 mm and a length of 
12.5 /zm for the interaction volume. The emitted spectral 
photon number is shown in Fig. [3j An integral over the 
spectrum yields an emitted photon number of 2.5 x 10~ 8 
at 50 keV per one collision of APTs which corresponds to 
a measureable signal of about 2 photons per day at 1 kHz 
repetition rate. Note that the choice of the atomic species 
is rather flexible. Multielectron highly charged ions offer 
an enhanced recombination probability due to core po- 
larization 23J but produce a larger electron background 
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FIG. 4: Schematic illustration of the chirping factor duin / dti ■ 
the harmonic bandwidth per unit ionization time. The insets 
show harmonic spectra for laser fields (black line) of either 
low (a) or high (b) intensity. The ionization time window 
At; resulting in HHG is marked on the time-axis, where the 
colors indicate the harmonic frequency originating from the 
particular ionization time. At; remains unchanged under the 
increase of intensity but the bandwidth of the contained har- 
monics increases from (a) to (b), consequently, decreasing the 
ionization propability per harmonic frequency expressed by 
an increase of the chirping factor. 



as duju/dti ~ wjj/Atj. Thus, the photon emission rate 
in a constant bandwidth for a zero-range potential scales 
as N n ~ Ip/tUtf 5 . A rough estimate for the scaling of I p 
at a constant ionization rate can be derived fixing the 
common tunneling exponent yielding I p ~ E 2 ^ 3 ~ <^f/ 3 
and, consequently, N„ ~ 1/wf^ 17 . The decrease for a 
hydrogen-like potential is even more dramatic. There- 
fore, the HHG photon yield decreases with rising photon 
energy due to the decreased probabilities of ionization 
per harmonic and the reduced recombination cross sec- 
tion. Our analysis points out a possible future direction 
for optimization of HHG by means of increasing the ion- 
ization time window at a given harmonic bandwidth. 

In conclusion, we have shown that the drift and phase- 
matching do not restrict HHG to the non-relativistic 
regime. The proposed setup renders the relativistic 
regime of HHG in a multi-atom ensemble accessible. 



that can be balanced by a lower gas density. The over- 
all efficiency is maintained or could even be enhanced. 
The bandwidth of phase-matched HHG in this scheme is 
about 150 eV near the cutoff and pulses with a duration 
of about 35 as can be produced. 

The small magnitude of the harmonic signal com- 
pared to current XUV HHG yields can be explained 
by investigating the spectral HHG photon rate N n for 
phase-matched emission [T3j [21] of the harmonic order 
n = u>h I uj from a fixed volume 



N n ~ Wi(ti) |(0|Vff|p)r Kr^WH/ffij) 



(5) 



Here ty$(i,-) is the ionization rate with the ionization time 
ti, (0|Vff|p) the recombination amplitude and the last 
factor accounts for the dynamical properties of the wave 
packet, Vj_t 2 expresses the transversal electron spreading 
with transversal spreading velocity v±, r the excursion 
time of the electron and du>H / dti is the so-called electron 
wave packet chirping factor discussed below. 

We proceed by analyzing the scaling of N n with in- 
creasing laser intensity at a harmonic energy near the re- 
spective cutoff provided that Wi(ti) is kept constant by an 
appropriate choice of I p . The recombination amplitude 
decreases with increased electron energy favoring scat- 
tering rather than recombination. Its scaling depends on 
the shape of the ionic potential: |(0|Vjj |p)c| 2 ~ Ip^ 2 I uj 'h 
for a hydrogen-like ion and |(0|Vff |p)z| 2 ~ ^/Lp/w 2 H for 
a zero-range potential with I p <C u>h and p 2 ~ ojh- Re- 
garding the last term of Eq. ([5]), we follow [3T] to find 

fx = VE/Ip^ 4 ~ (ojh I Ip) 1 ^ and illustrate the chirping 
factor in Fig. |4j It describes that the bandwidth of the 
harmonics emitted from a fixed ionization time window 
rises with increasing laser intensity (i.e. the ionization 
probability per harmonic decreases) and can be estimated 
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